INLO-PUB-18/94 

The sphaleron energy for SU(2)-Higgs from cooling 

Margarita Garcia Perez^ and Pierre van Baal ^ * 

''Instituut-Lorentz for Theoretical Physics, 
University of Leiden, PO Box 9506, 
NL-2300 RA Leiden, The Netherlands. 

The cooling algorithm for saddle points presented in ref. [1] is generalized to obtain static classical solutions of 
the SU(2)-Higgs field theory in the limit of infinite Higgs self-coupling. The sphaleron energy obtained via this 
algorithm is £'spi, = 5.08(7) Mw /ctw to be compared with 5.0707 obtained in the variational approach [2]. 
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The spherically symmetric sphaleron solutions 
for the SU(2)-Higgs field theory, obtained by a 
variational analysis, have been known for quite 
some time [3,4]. Above Mh = 12Mw the spha- 
leron undergoes a series of bifurcations [2], ac- 
quiring at each bifurcation an additional neg- 
ative mode, while a new solution, a so-called 
deformed sphaleron splits off. For infinite A, 
where the model is identical to the gauged non- 
linear sigma model, there is an infinite number 
of solutions ranging in energy from SAlMw/ctw 
to the energy of the lowest deformed sphaleron 
5.07 Mw/ ctw 7 which has only one negative mode 
with Lu^ = — 4.714M^ (the number of unstable 
modes increases with increasing energy). At in- 
finite Higgs self-coupling these solutions are re- 
lated to the so-called electro- weak skyrmions [5]. 
Generalizing the pure gauge cooling algorithm for 
saddle points, we present the lattice results for the 
energy of this lowest deformed sphaleron, hence- 
forth called the sphaleron (using the more restric- 
tive definition that requires the existence of pre- 
cisely one unstable mode). We view it as a useful 
check for the efficiency of the algorithm, but one 
might envisage useful applications for the compu- 
tation of sphaleron transition rates [6]. 

1. The algorithm 

The dynamical variables for the SU(2)-Higgs 
model on the lattice are the gauge group vari- 
ables Vfi{x), defined on the link that runs fron x 



to X + fi, and the Higgs field in the fundamental 
representation of SU(2). For infinite self-coupling 
in the Higgs sector, the length of the Higgs dou- 
blet is frozen and can be chosen equal to unity. 
The Higgs field can in this case be represented by 
a SU(2) matrix crix), which is associated to the 
gauge degree of freedom and can be reabsorbed 
into the links via the change of variables [7] 

Uf,{x) = a{x)Vf,{x)a{x + ij) . (1) 

The lattice action is (C7^(a;) = •-^) 

*=|?|Et.(i--D)+Ct.(i-1)|-(2) 

For n = 4 (n = 3) the continuum action (energy) 
functional is recovered by setting k = 2M^a^, 
with a the lattice spacing. In what follows we 
will restrict to the n = 3 case, but the method is 
applicable in any number of dimensions. 

Designing a cooling algorithm for this model 
goes exactly as for the pure gauge theory. We 
construct the functional S by squaring the equa- 
tions of motion for S 

S=4^Y.'^^(u,{x)Ul{x)-{U,{x)Ul{x)f) , (3) 



X,IJ, 



where 



^.w=?+E('^n+o 



(4) 



vi^ij. 



* Based on the talk presented by the first author at Lat'94 
(Bielefeld, Sept. 1994) 



The equations of motion for S are solved by 
UM = MJ^{x)l\\lJ,{x)\\ . (5) 



The + sign is to be taken in order for the solution 
to have a smooth continuum limit. Minimizing S 
amounts to iteratively replacing C/^(a;) by eq. (5). 
To minimize S with respect to a single link vari- 
able would require us to solve an eighth order 
polynomial in ||C7p(a;)|| at each iteration. An al- 
ternative algorithm [1] can be designed by noting 
that S will always be lowered under the update 



U'{^) = 



M{U^{x))-W^{x) 



\\M{U,{x))-W,{x)\\ ' 
with {V"{x) are the 2(n — 1) staples in eq. (4)) 

M{U,{x)) ^ 2 Tr (u,{x)Ul{x)) U,{x) 
+&Y.^v{U,(x)V^(x)^)V^(x) 



(6) 



(7) 



and 



i'w=iEci-:i)n+E!-- 
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with the unit vector a G {±1, ■ ■ ■ , ±n}. We give 
the explicit form for the term proportional to k, 
refering for the other terms to ref. [1], 



nY,iI-U!{x))U,{x + a)UUx + fi) 



(8) 



with the convention U-a{x) = U^{x — a). This 
algorithm will have any extremum of the energy 
functional, even lattice dislocations, as a mini- 
mum. To avoid being trapped in one of these 
dislocations we first bring S down with eq. (5) 
to a smooth configuration. We then use eq. (6). 
When this does not lower S further some updates 
with eq. (5) often help (see ref. [1] for details). 

2. Results 

Before presenting the results it is useful to clar- 
ify a few points. We are looking for infinite vol- 
ume solutions. To minimize boundary effects due 
to the finite lattice size, the typical correlation 
length of the system has to be much smaller than 



the size of the box (aMw)^^ = \/2/k <ti N, on 
a lattice of size N^ . This gives, for fixed N, a 
lower bound on the values of k to be considered. 
The simulations were performed for N = 8, 12 
and 16, imposing A^a/k/2 > 2.5. For smaller k 
values it turns out that the electro-weak sphale- 
ron developes additional unstable modes. This is 
not due to a bifurcation, but due to the crossing 
in energy of two widely separated critical points. 
The other critical point responsible for this is the 
finite volume sphaleron, constructed for k = 
in ref. [1]. (It is interesting to note that this is 
not forbidden by Morse theory, which measures 
the global topology. Only if the Morse functional 
- i.e. the energy functional - has no degenerate 
critical points, bifurcation is implied when addi- 
tional unstable modes appear). 

On the other hand, for large k at fixed A'^, 
the energy density becomes highly peaked over a 
few lattice spacings, giving large lattice artefacts. 
One would like the correlation length to be much 
bigger than the lattice spacing, a/2/k ^ 1. In the 
continuum there is an infinite set of sphaleron-like 
solutions in a very small range of energies, with 
no bound on the number of unstable modes. On 
the lattice there can, however, not be more un- 
stable modes than the finite number of degrees 
of freedom. The way the solutions bifurcate now 
also depends on the lattice spacing. The largest 
lattice spacing we can allow is the one where our 
sphaleron acquires more than one unstable mode. 
This can be monitored by computing the Hessian, 
which is only practical for N = 8. From this we 
determined the bound that N ^J njl < 3.2 (i.e. 
the correlation length has to be slightly bigger 
than two lattice spacings). 

The way we obtained the required configura- 
tions was by starting at A'^ = 8 with the links at 
the boundary frozen to unity. This lifts the en- 
ergy of the finite volume sphaleron considerably, 
such that the value of k at this stage may even 
be below the value we quoted above. Also, lattice 
artefacts cause the breakdown of translational in- 
variance on a periodic lattice, thereby generating 
spurious saddle points with up to three extra un- 
stable modes. It has the additional advantage of 
centering the energy profile in a maximally sym- 
metric way, whereby one minimizes the lattice 




Figure 1. fgph as a function of a'^M^ fitted for 
fixed physical volume to a linear function of a^. 

artefacts. We then release the frozen boundary 
condition and compute after cooling the Hessian 
to verify that we have one unstable mode only (at 
this stage k has to be brought to the window of 
values discussed above). For A^ = 12 and 16 the 
initial configurations were generated from the one 
at A^ = 8, by embedding it in the large lattice 
(links parallel to the boundary remain constant 
and those perpendicular to the boundary are put 
to unity) and adjusting k such that the physical 
volume is at least as large as for A'^ = 8. The 
results for the sphaleron energies, as well as the 
negative eigenvalue for the N = 8 Hessian are 
presented in the table below. The values of S are 
not larger than 10~^M^/aw and even consider- 
ably smaller for the smaller lattices. 




Figure 2. Continuum extrapolated values for fgph 
as function of the physical volume MwL. 

the three pairs of lattices) are plotted as a func- 
tion of the physical size. Assuming that one can 
compute the fields in terms of the free massive 
propagator outside of the core region, one easily 
derives the volume dependence to be 

4ph(i) = 4^h + 7 6-^-^(1 + OiMwL)-') . (9) 

Using the extrapolations linear in a^ we find 
for E'^-^ a value of 5.02(l)Mw/aw, dropping one 
(on which the fit in fig. 2 is based), two or three of 
the smallest volumes gives resp. 5.04(1), 5.05(1) 
and 5.06. On the other hand, using an extrapola- 
tion quadratic in a^ gives resp. 5.01(3), 5.08(3), 
5.12(3) and 5.16. As a conservative estimate we 
give £^^ = 5.08(7) Mw /ctw^ covering all values. 
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To get rid of lattice spacing errors we first work 
at fixed physical volume set by MwL = N \J nf^ 
for three values of the lattice spacing a = L/N. 
We extrapolate a — > by fitting the data to a lin- 
ear function of a^ (fig. 1) . To check the stability of 
the fit we also performed a fit quadratic in a^. In 
fig. 2 the linearly extrapolated values (errors de- 
termined by averaging intercepts computed from 
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